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Abstract 

This is the second in a series of two lectures on the technique of dimensional continuation, a new 
method for analytically calculating certain energy transport quantities in a weakly to moderately 
coupled plasma. Recently, this method was employed by Brown, Preston, and Singleton (BPS) to 
calculate the electron-ion temperature equilibration rate and the charged particle stopping power 
to leading and next-to-leading order in the plasma coupling. The basic idea is very simple. Con- 
centrating upon the equilibration rate, the calculation consists of the following two steps: (i) per- 
turbatively expand the rate in the form d£/dt = —Ag 2 In g + Bg 2 + 0(g 3 ), with the dimensionless 
expansion parameter being defined by g = e 2 K e /4irT e ; (ii) analytically calculate the coefficients A 
and B using the method of dimensional continuation. The factor of 4-7T should be omitted from g in 
nonrationalized electrostatic units. In the first lecture, I presented a basic overview of the requisite 
theoretical machinery of dimensional continuation imported from particle physics, but in a self- 
contained manner that assumed no familiarity with quantum field theory. In this lecture, I develop 
the framework further, and then explicitly calculate the electron-ion temperature equilibration rate 
in the high temperature limit. In this extreme quantum limit, the calculation of the coefficients A 
and B simplifies considerably, allowing us to concentrate on the physics of the method rather than 
the added complexity of the more general BPS calculation. This method captures all short and 
long distance physics to second order in g, while three-body and higher correlations are contained 
in the cubic and higher order terms denoted by 0(g 3 ). In a weakly to moderately coupled plasma, 
where g is small, the error term 0(g 3 ) in this calculation is also small compared to the A- and 
i3-terms, in which case the BPS methodology is quite accurate. Should higher order contributions 
be required, they can be calculated systematically, thereby improving the accuracy of the result in 
a controlled manner. To get a feel for the numbers, one finds g ~ 0.04 at the center of the sun, 
where the plasma conditions are n ~ 5 x 10 25 cm -3 and T ~ 1 keV. The coupling constant g can be 
scaled to other plasma regimes through the proportionality relation g oc n 

1/2 2-3/2. 0f 

course the 

application of interest determines the relevant plasma regime, which may or may not lie within the 
domain of applicability of the BPS calculation. For example, the technique breaks down for warm 
dense matter where is g not very small; however, this analytic perturbative technique is applicable 
for ignition in inertial confinement fusion and for other processes in hot a weakly coupled plasma. 
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I. INTRODUCTION AND REVIEW 



This is the second lecture on dimensional continuation, a new technique jl] recently used 
to calculate the charged particle stopping power and the temperature equilibration rate in a 
weakly to moderately coupled plasma [21]. In Lecture I [3j of this series, I discussed the basic 
theoretical machinery of dimensional continuation, and I performed a model calculation of 
the equilibration rate. Reference [4] also contains a summary of the method in a very readable 
form. In this lecture, I will present the complete calculation of the electron-ion temperature 
equilibration rate in the extreme quantum limit, valid to leading and next-to-leading order 
in the number density (a more general calculation is performed in Section 12 of Ref. [2j to all 
orders in quantum mechanics, thereby providing an exact interpolation between the extreme 
classical and quantum limits). This calculation is near exact for a weakly coupled plasma, 
and it is quite accurate for a moderately coupled plasma. Before proceeding directly to the 
calculation, however, it might be useful to quickly review some of the more salient features 
of dimensional continuation discussed in Lecture I. 

Under most circumstances, a plasma is not produced in thermal equilibrium; for example, 
when a laser ionizes a substance, it preferentially heats the electrons over the ions. However, 
since the electrons are so light, they rapidly come into thermal equilibrium among themselves 
with temperature T e ; some time later, the ions too will equilibrate among themselves to a 
common temperature T : . Finally, the electrons and ions will begin to equilibrate, and it is this 
process upon which we shall focus. Let d£ eI /dt denote the rate per unit volume at which the 
electron system at temperature T e exchanges energy with the ion system at temperature T z 
through Coulomb interactions (throughout these notes, I will always measure temperature 
in energy units). The electron- ion equilibration rate is proportional to the temperature 
difference between the electrons and ions, and can be expressed by 

^i = -C eI (T e -T I ) . (1.1) 

To restate the goal of this lecture more precisely, we shall calculate C eI in the high temperature 
limit [where two-body scattering is accurately given by the Born approximation], and we 
will do so exactly to leading and next-to-leading order in the plasma coupling parameter g 



[defined in Lecture I, or in 
a particularly simple form 



q. (12. II) of this lecture]. Under these conditions, the result takes 



where 7 = 0.57721 • • • is the Euler constant, n e and u e are the electron Debye wave number 
and plasma frequency, and ojf = J2i u i * s sum °f the squares of the ion plasma frequencies. 1 
In the form displayed by equation fll.2p . the rate coefficient C el and the Coulomb log- 
arithm lnA BPS do not explicitly depend upon one's choice of electrostatic units, and one 
may calculate the Debye wave numbers and the plasma frequencies in any desired system. 
For dimensional continuation, however, it is more convenient to use rationalized electrostatic 
units, and I shall employ this choice from here out. An arbitrary plasma component will be 
labeled by an index b, and is characterized by mass m b , charge = Z b e, number density 
Ub, and temperature TJ,. The index b can span the electron and ion plasma components, 
that is to say, b — e,i with i being an arbitrary ion species. Working in three dimensions 
for now, the Coulomb potential between two charges e a and e& separated by a distance r is 
V = e a eb/4:7ir, and in rationalized units, the Debye wave number and the plasma frequency 
of species b take the form 2 

= ^ (1-3) 

ut= e J^L. (1.4) 
The square of the total Debye wave number is k 2 = J2b K &> an d the total Debye wave length 



is A D = k d 
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II. CALCULATING THE RATE IN PERTURBATION THEORY 



Reference [2[, hereafter referred to as BPS, used a double pronged strategy to calculate 
the rate coefficient (11.21) . First, a well chosen [5J dimensionless parameter g was constructed 
from the relevant dimensionfull plasma quantities, thereby providing a parameter in which to 
perform a controlled perturbative expansion to leading and next-to- leading order in g. The 
systematic error of the calculation was estimated by the cubic order term in the expansion, 
which is quite small for a weakly to moderately coupled plasma. While perturbative calcu- 
lations are not very common in plasma physics, primarily because of the complexity of the 
systems of interest and the computational focus within the field, the validity of perturbation 
theory should nonetheless be clear for a "simple" system such as a weakly coupled and fully 
ionized plasma. The second part of the BPS argument deployed a powerful technique from 
quantum field theory allowing one to analytically calculate the coefficients in the ^-expansion. 



1 Equation (|1.2p corresponds to Eqs. (3.61) and (12.12) of Ref. [2j, where I have taken this opportunity to 
correct a small transcription error: when passing from Eq. (12.43) to Eq. (12.44) in Ref. [2|, a factor of 
1/2 was dropped. Restoring this factor of 1/2 changes the additive constant outside the logarithm from 
the —7 — 2 that appears in Eq. (12.12) of Ref. Q to the constant —7 — 1 in (|1.2| . 

2 In nonrationalized units, the right-hand-sides of (|1.3p and (|1.4|) should contain an additional factor of in. 
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A. Perturbative Expansions in Weakly Coupled Plasmas 

Let us first concentrate on the perturbative expansion. As demonstrated in Ref. |||, and 
discussed at length in Lecture I [3j, for the case at hand the dimensionless plasma coupling 
parameter is defined by 3 



Note that g is the ratio of the Coulomb potential energy of two point-charges, separated by 
the screening length k" 1 , to the thermal energy of the plasma. Therefore, g can be used 
to measure the strength of the plasma. To get a feel for the size of g in a hot but not too 
dense plasma, one finds g = 0.042 for a hydrogen plasma under the solar-like conditions 
n e = 5.0 x 10 25 cm -3 and T e = 1.3 keV. One can scale to other density and temperature 
regimes by noting that g oc n e /2 Te 3/2 . It was shown in Ref. [s| that plasma quantities 
always expand in integer powers of the coupling g, and therefore g is the appropriate pa- 
rameter in which to perform a controlled perturbative analysis for weakly coupled plasmas. 4 
The ^-expansion allows for possible non-analytic terms, such as In g, and in particular, the 
electron-ion equilibration rate can be written 



where I have indicated the leading order (LO) and the next-to-leading order (NLO) terms in 
the expansion. The minus sign on the leading order term of (12.21) is a matter of convention, 
and for small values of g it renders the coefficient A positive when the energy exchange is 
positive. Provided we can calculate the coefficients A and B, then (12.21) will be quite accurate 
in a weakly to moderately coupled plasma in which g is small. Of course this perturbative 
approach breaks down for strongly coupled plasmas, those for which the value of g is of order 
one or greater, since every term in the expansion becomes equally important in such cases. 
However, unlike a model or an uncontrolled calculation, the BPS calculation informs us of 
its domain of validity, and it provides an estimate of its own error through the size of g. 

3 In nonrationalizcd units we would write g = e 2 K e /T e , with k 2 = 4ir e 2 n e /T e . 

4 The usual plasma parameter T is related to the expansion parameter by g oc r 3 / 2 (with proportionality 
constant of order unity). Small values of T therefore imply small values of <?, and one may characterize 
the strength of the plasma by either g or T. The proportionality relation above follows from the fact 
that g oc n X J 2 and V oc n^ 3 (the parameter g is defined in terms of a Debye screening length rej 1 , while 
r is defined in terms of the inter-particle spacing n e 1 ^ 3 ). Furthermore, since g oc nj" 1 , we may loosely 
think of the ^-expansion as an expansion in the electron number density, as I have done in the first 
paragraph of this introduction. More precisely, of course, we are expanding in the dimensionless quantity 
g oc e 3 n 

i/2 T - 3 /2 gee Ref _ g for 

more details, particularly Section 1.1 entitled Relevant Scales and 

Dimensionless Parameters. 




(2.2) 



LO 



NLO 
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B. Calculating the Coefficients of the Expansion 



We have now reduced the problem to finding the coefficients A and B of the rate (12. 2p . 
The coefficient A was first obtained long ago by Spitzer (and it can be estimated by 
dimensional analysis alone). The coefficient B, however, was calculated only recently in 
Ref. [2J, which employed a powerful technique from quantum field theory called dimensional 
regularization, or dimensional continuation as I will call it here. Since this technique is quite 
subtle and has proven to be somewhat controversial, I should emphasize that the method by 
which one chooses to calculate these coefficients is immaterial, except to the extent that it 
must contain enough physics to extract the next-to- leading order coefficient B. Techniques 
other than dimensional continuation could well furnish one with the correct expressions 
for A and B, and perhaps in a simpler manner. However, the only relevant point here is 
that by hook or by crook we must analytically calculate these coefficients, and dimensional 
continuation is one method of doing this. 5 

Before turning to the calculation of the coefficients, allow me to make a comment on 
the relation between the next-to-leading order B-term and the Coulomb logarithm. Writing 
the leading order coefficient as K = Ag 2 , and defining the dimensionless coefficient C = 
exp{— B/A}, we can express the rate (I2.2p in the form 



Since the Coulomb logarithm means different things to different people, 6 I would like to be 
quite specific in this lecture. By the words "Coulomb logarithm" I simply mean the term 
lnA coul defined in ( 12.31) . excluding the cubic and higher order terms. Hence, calculating the 
next-to-leading order coefficient B is equivalent to determining the dimensionless coefficient 
C inside the Coulomb logarithm. Finding dimensionless constants is usually a difficult 
problem, particularly since one cannot appeal to dimensional analysis for an estimate. It 
should therefore not be surprising that the coefficient C varies over an order of magnitude 
or so across the various models within the literature. 



5 I have recently been informed Q that Ref. [7( , which is designed to apply to both strongly and weakly 
coupled plasmas, reproduces the BPS result (jl.2|) in the limit of weak coupling. As far as I am aware, 
Refs. @ and 0] are the only works currently in the literature with a formalism strong enough to extract 
such next-to-leading order physics from first principles. 

6 Student: What is the Coulomb logarithm? Professor: 10. 




with In A, 



coul 



In {CM 



(2.3) 
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III. CALCULATING IN ARBITRARY DIMENSIONS 



Before proceeding directly to the calculation in Sec. HVl let us further develop the basic 
physics and mathematical machinery necessary to perform calculations in an arbitrary num- 
ber of dimensions. The motivation for this section is, of course, a thorough exposition of 
the BPS methodology for calculating Coulomb energy-loss processes in a plasma. However, 
the material in this section is well known and applicable to a wide variety of other calcu- 
lations, such as particle decay rates in high energy physics and analytic work in statistical 
mechanics. For the sake of completeness, however, and to establish some results that will be 
useful in Sec. HVt I will present a cursory but self-contained treatment here. If this material 
is familiar, then one may proceed directly to the calculation of the temperature equilibration 
rate in Sec. IIVI (given the background material in this section, the calculation itself is less 
than eight pages in length). 

We shall start by developing the hyperspherical coordinate system in v dimensions, which 
is a straightforward generalization of three dimensional spherical coordinates. To illustrate 
the utility of hyperspherical coordinates, I will calculate the hyperarea and hypervolume of 
several multidimensional objects by exploiting their spherical and cylindrical symmetries. 
These results will be used quite extensively in the next section. As a physical application, 
I then develop the multidimensional analog of the scattering cross section, which will al- 
low us to consistently include short- distance quantum scattering effects in the ^-expansion 
(quantum effects manifest themselves through the ^-dependence of the coefficients in this 
expansion). Since we are interested in Coulomb energy exchange, we next examine electro- 
statics in arbitrary dimensions. From the multidimensional form of Gauss' Law, we shall 
derive the //-dimensional Coulomb potential Vj,(x), and we will see that it depends only 
upon r = |x| in such a way as to emphasizes short distance physics when v > 3 and long 
distance physics when v < 3. In v = 3, the short and long distance physics compete with 
equal strength, giving an infrared and an ultraviolet divergence, and this is what renders 
the temperature equilibration problem so difficult. To employ the extreme quantum limit, 
in which the Born approximation for the two-body scattering dominates, we must calculate 
the Fourier transform of the Coulomb potential in v dimensions. Interestingly, we shall find 
that the Fourier transform of V u (r) is given by the quite simple expression K,(k) = 1/k 2 , the 
form of which does not depend upon the dimension of space, but only upon the length of the 
wavenumber k — |k|. The fact that V v {k) is so simple greatly facilitates calculations in the 
extreme quantum limit. With potential in hand, we shall then construct kinetic equations 
in v dimensions. These equations are explicitly finite in all but v — 3 dimensions, and I will 
explain the manner by which the BBGKY hierarchy reduces to the Boltzmann equation and 
the Lenard-Balescu equation (in v > 3 and v < 3 respectively). 
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A. Kinematics and Hyperspherical Coordinates 



1. Hyperspherical Coordinates 

Kinematic quantities such as the //-dimensional momentum or position vectors are ele- 
ments of the same //-dimensional Euclidean space W. For definiteness, I will specialize to 
the case of position r, with the understanding that this vector could also refer to momentum 
or wavenumber. We can decompose any vector r e W in terms of a rectilinear orthonormal 
basis e£, so that r = Y^e=i x e ^ or m component notation r = (xi, ■ ■ • , x v ). Each component 
is given by xg = eg ■ r, and a change dr in the vector r corresponds to a change dxg = eg • dr 
in the rectilinear coordinate xg. Letting r vary successively along the independent directions 
eg, we can trace out a small //-dimensional hypercube with sides of length dxg\ therefore, the 
rectilinear volume element is given by the simple form 



In performing integrals over the kinematic variables, however, symmetry usually dictates 
the use of hyperspherical coordinates rather than rectilinear coordinates. I will therefore 
review the hyperspherical coordinate system in this subsection, deriving the measure for 
a //-dimensional volume element d v x in terms of hyperspherical coordinates. For our pur- 
poses, the primary utility of hyperspherical coordinates is that the volume element d u x can 
be written as a product of certain conveniently chosen dimensionless angles, which I will 
collectively refer to as dQ v _i, and an overall dimensionfull radial factor r u ~ 1 dr, so that 
d v x = dn v _ ir u ~ l dr. 

Starting with the usual 3-dimensional spherical coordinates of Fig. [TJ let us recall why 
the three dimensional volume element takes the form d 3 x = sin 9 dd dip r 2 dr (with < 9 < it 
and < < 2tt, and of course < r < oo; the coordinate singularities of the spherical 
system are not important here). As depicted in the figure, the three dimensional vector r 
has length r, and subtends a polar angle 9 relative to the z-axis, while its projection onto 
the x-y plane subtends an azimuthal angle relative to the x-axis. The two angles 9 and 
specify completely the direction of the unit vector f . As we increase the polar angle 9 by a 
small amount dO, the vector r sweeps out an arc of length dR\ = rd9; similarly, a change dcf) 
in the azimuthal angle will cause r to sweep out a perpendicular arc (in the x-y plane) of 
length dB,2 = rsm9dcf). Note that the factor of sin^ in dR,2 arises from the projection of r 
onto the x-y plane. We can make one more independent displacement by moving dr units in 
the radial direction, which results in a line of length dR 3 = dr. For small displacements in d9, 



u 



d v x = J_ J_ dxg = dx\ dx2 • ■ ■ dx, 
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X 

FIG. 1: Spherical coordinates r, 9, <f> of a point r in three dimensional space: radial distance r, polar angle 
9, and azimuthal angle <f>. The angles range over the values < 9 < ir and < 4> < 27r. 




(a) (b) 

FIG. 2: Hypcrspherical coordinates r, #1,6*2, #3 of a point r in four dimensional space. As before, r = |r| is 
the radial distance. The angles are defined as follows, (a) First, let 9\ be the angle between r and the w-axis. 
Let us now project r onto the orthogonal three dimensional space, so that r = [x, y, z, w) — > r w = (x, y, z, 0). 
The length of this projection is r w — rsin^i, and the projection itself is the same as projecting r onto 
the three dimensional hyperplanc w — 0. (b) The vector r w can be viewed as a three dimensional vector 
r w = (x,y, z), which then defines the usual polar and azimuthal angles of Fig.[TJ denoted here by #2 and #3 
respectively. 

d<p, and dr, the vector r sweeps out a small cubic volume element with sides of length dRi, 
dR<ii and dR%. The volume of this element is therefore d 3 x = dRi dR<i dR% = rd9-r sin Odip-dr. 

Let us now consider the volume element d A x in four dimensional space. Denote the 
coordinates of a vector r by x,y,z,w, that is to say, take r = (x,y,z,w). Since we cannot 
visualize four dimensional space, 7 let us examine this problem in two steps, each of which 

7 Apart from visualization problems, we can nonetheless work in higher dimensions by employing analytic 
geometry and analogies with lower dimensions. For example, an ordinary two-sphere of radius r, which 
I will denote by S^fV); has the equation x 1 + y 2 + z 2 = r 2 in three dimensional space; a corresponding 
"three-sphere" S%(r) in four dimensions can be represented by x 2 +y 2 + z 2 +w 2 = r 2 . In a similar manner, 
a two-dimensional plane in three-space can be expressed as a\ x + a.2 y + 03 z — c for real numbers ae 
and c, while a three-plane in four dimensional space takes the form a\X + azy + a$z + d^w = c. As 
a final example, consider a "three-cone" oriented along the w-axis: x 2 + y 2 + z 2 — w 2 = 0. The "conic 
sections" are obtained by slicing this hypercone with a three-plane along various orientations. For example, 
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can be visualized in either two or three dimensions. First, consider the plane that contains 
the w-axis and the vector r, and let 9i be the angle between the u>-axis and the vector 
r in this plane, as shown in Fig. [2^. We now project r onto the w = hyperplane (a 
three dimensional slice of four-space), calling the projected vector r w . Since the three-plane 
w = lies perpendicular to each of the axes x, y, and z, the vector lies in the three 
dimensional space shown in Fig. [2b, and its length is \r w \ = rsin^x. Let the angle 82 be the 
polar angle between the z-axis and the vector r w , while #3 is the usual azimuthal angle 0, 
as illustrated in Fig. [2b. As we vary the three angles and the radial coordinate, we sweep 
out a four- dimensional cube (or an approximate cube) with sides of length dR\ = r d9±, 
dR,2 = r sin#i<i#2, dR% = r sin 9\ sin 62 d9 3 , and dR A = dr. This gives a four dimensional 
volume element 

d 4 x = dR x dR 2 dR 3 dR A = sin 2 9 x d9 x sin 9 2 d9 2 d9 3 r 3 dr , (3.2) 

where < 0£ < ir for £ = 1,2 and < #3 < 2n. As a useful exercise, we can find the four 
dimensional hypervolume enclosed by a three-sphere of radius r by integrating the volume 
element over the appropriate bounds, 

B A = I d9 l sm 2 9 l j d9 2 sm9 2 / d9 3 / dr'r' 3 = - vr 2 r 4 . (3.3) 
Jo Jo Jo Jo 2 

The derivative of B4 with respect to r gives the hypersurface area of the three-sphere, 

S 3 = d J± = 2^r>. (3.4) 

This is analogous to a three dimensional ball of radius r and volume B 3 = 47rr 3 /3 bounded 
by the two-sphere of area S2 = 47rr 2 . 

We can readily generalize this procedure to an arbitrary number of dimensions. Consider 
a point rel 17 given by the rectilinear coordinates r = (xi, x 2 , ■ ■ ■ , x v \ Let 9\ be the angle 
between the vector r and the Xi-axis, in a manner similar to that of Figs. [1] and [2^. Note that 
dR\ = rd9\ is the arc length swept out by r as the angle 9\ is incremented by dB\. Let us now 
project r onto the hyperplane x\ = 0, the {y — l)-plane normal to the Xi-axis and passing 
through the origin, calling this projection r x : that is to say, let r — > r 1 = (0,x 2 , • • • ,x u ). 
The length of this vector is r\ = rsin^. Let us proceed to the next step and define the 
angle 9 2 as the angle between the avaxis and r 1; in which the angle 02 is varied 

by d9 2 , the vector r x sweeps out an arc of length dR 2 = rid9 2 = rsm9id9 2 . In a similar 

if we slice the three-cone by the hyperplane w — r orthogonal to the w-axis, then we find a two-sphere 
x 2 + y 2 + z 2 = r 2 ; if we slice the three-cone by a hyperplane along the z-axis, say z = r, then we find the 
hyperboloid of two sheets x 2 + y 2 = w 2 — r 2 oriented along the w-axis. 
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TABLE I: Solid angle f2„_i as a function of dimension v. 



V 


2 


3 


4 


5 


6 


7 


8 


20 




2tt 


4tt 


2^ 2 


8vr 2 /3 


7T 3 


16vr 3 /15 


7T 4 /3 • 


vr 10 /181440 


value 


6.28 


12.6 


19.7 


26.3 


31.0 


33.1 


32.5 


0.516 



fashion, project ri onto the X2-plane, that is, the plane described by x\ = and %2 = 0. 
This projection is given by r — ► r 2 = (0,0, £3, • ■ ■ ,x u ), and the length of the projection is 
T2 = rising = t sin 9\ sin Qi- 8 For the general £ th iteration, let 8e be the angle between 
the X£-axis and r^_i, so that dRe = re~\ d9i = r sin 9\ sin 82 ■ ■ ■ sin 6^_i d^. In summary, we 
define the quantities 

6e = angle between the x^-axis and r^_i (3.5) 
dRe = re_i d9% = r sin 9\ sin 62 ■ ■ ■ sin Bg-\ dBg (3-6) 
r -> r t = (0, • • • ,0,x e+1 , • • • ,x„) (3.7) 
r£ = ri_i sin Qi = r sin #1 sin 62 ■ ■ ■ sin Qi-\ sin 6^ , (3-8) 

where we have used the fact that re-i = r sin d\ sin #2 • ■ • sin 9e-i- The last two are lines 
provide the projection for the (£ + l) st step. This gives the //-dimensional volume element 

V 

d v x = Y[dRt = sin"- 2 B x dQ x ■ sin 1 '" 3 6 2 d0 2 ■ • • sin B v - 2 d6 v - 2 ■ dQ v - X ■ r^dr . (3.9) 

For notational convenience, I will write the angular measure in (13.91) as dO„_i, so that 

d v x = dn v -ir"- l dr . (3.10) 
As we proved in the Lecture I, the integration over all angles gives the total solid angle 

r 2W 2 

<Vi = J = —j^r , (3.11) 

and Tabled illustrates the numerical values of this solid angle over a wide range of dimensions. 
Note that fl u -i reaches a maximum around v — 7 and then slowly decreases. 

8 Just for good measure, let us do one more iteration. Let 83 be the angle between the a^-axis and the 
previous projection T2, in which case (IR3 = ^(103 = rsin#i sin 82 (I83. Let us now project T2 onto the 
a;3-plane described by x\ = 0, X2 = 0, and X3 = 0, i.e. r — > r3 = (0, 0, 0, £4, • • • , x v ). The length of this 
vector is r3 = r2 sin #3 = r sin#i sin #2 sin $3, and the next iteration can begin. 
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As a matter of completeness, let us prove (13.111) here. First, consider the one- dimensional 
Gaussian integral 



oo 

2 



die" 1 = v 7 ^ • (3.12) 
If we multiply both sides together v times (with u £ Z + ), we find 

/CO pOO p OO /* 

dx l e~ x * dx 2 e- x '--- dx u e~ xl = d u xe~ r \ (3.13) 
-oo J —oo J —oo J 

where the vector r in the exponential of the last expression is the i/-dimensional vector 
r = (xx,x 2 , ■ ■ ■ ,x u ), and r 2 = r ■ r = Y^i=i x £ ■ As in (13-101) . we can factor the angular 
integrals out of the right-hand-side of (I3.13p . and the remaining one-dimensional integral 
can be converted to a Gamma function with the change of variables t = r 2 : 

ir u/2 =Jdn u -t -J drr v ~ l e~ r2 =Jdtt v ^ ■ ^T(u/2) . (3.14) 

Solving for f dQ u ^ in (l3~T4l) gives d3~TT]) . 

In calculating the temperature equilibration rate and the charged particle stopping power, 
we encounter integrals of the form 



oo 

u-1 



h{v) = Jd v xh{r) = drr^ 1 f^r) (3.15) 

//•OO P7T 
d u xf 2 (r,6)= Q u _ 2 drr- 1 d6 sin- 2 6 f 2 (r, 6), (3.16) 
Jo Jo 

respectively, with v G Z + . The exact forms of f\ and f 2 do not concern us here, except that 
their angular dependence is determined by the following considerations: the integral (13.151) 
is spherically symmetric since the energy exchange between plasma species is isotropic, while 
in the latter integral (13.161) . the motion of the charged particle defines a preferred direction 
around which one must integrate. The integrals ii and I 2 can be viewed as functions defined 
on the positive integers, and as discussed at length in Lecture I p], Carlson's Theorem 
ensures that there is a unique analytic continuation onto the complex plane. As our first 
application in this section, let us see how the expressions (13.151) and (13.161) provide a means 
by which to easily and conveniently perform this analytic continuation to complex values 
of the spatial dimension u, thereby rendering v truly arbitrary. First, the solid angles fl u -i 
and il u - 2 are composed of a simple exponential factor tt u ^ 2 and a Gamma function, whose 
analytic continuations have been well studied. As for the integrals, simply treat v as a 
complex parameter, performing the one dimensional integral (13 . 15[) and the double integral 
(I3.16P in the usual manner of ordinary calculus. This provides functions I\{y) and l 2 {y) of a 
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By(r) 




FIG. 3: A (y — l)-dimensional sphere S u -\ of radius r bounds the z/-dimensional ball B u (r) of radius r. By 
integrating over successive shells of area, we can find the volume by B u {r) = J dr' Su—i(r'); or conversely 
S u (r) = B' v (r). 

complex argument v G C, in fulfillment of Carlson's Theorem. Double integrals of the form 
( I3.16P were used extensively in Ref. [2| to calculate the stopping power, where the angle 9 is 
determined by the direction of motion of the charged particle. Calculating the temperature 
equilibration rate, on the other hand, requires only the simpler one dimensional integral 
(I3.15p . as the energy exchange in this process is isotropic. 



2. The Hypervolume of Spheres, Disks, and Cylinders 

We shall now calculate the hypervolume of several useful geometric objects. Let us first 
consider a v- dimensional ball of radius r, defined by the set of points x G ~W for which 
|x| < r. We will denote this object by B u (r), and in two and three dimensions this is a 
disk and a spherical, both volume centered at the origin. We can find the v- dimensional 
hypervolume of the ball B v (r) by simply integrating (13.91) over all permissible values of the 
coordinates. It should cause no confusion to denote the hypervolume of the region B u (r) by 
the same symbol, and using (13.111) we find 

f r W 2 

B v {r) = dtt^ dr'r' v - 1 = — — r v . (3.17) 



J Jo I>/2 + l, 

The boundary of B u (r) is a [y — l)-dimensional sphere S u _i(r) defined by |x| = r, or 
Y^e=i = r 2 . By differentiating (13 . 1 T[) with respect to the radius r, we can also find the 
hyperarea of a (y— l)-dimensional sphere S u -i(r) of radius r in M. p , 



ii 



FIG. 4: The hyperarea of a hypercylinder C„_i(r, L) of length L and radius r is C„_i(r, L) = 5^-2 (?*) ■ £ f 
and the hypervolume bounded by the cylinder is V v {r, L) = B v -i{r) ■ L. 

For brevity, I have denoted the hyperarea by the same symbol S„-i(r) as the sphere itself, 
which is simply the (V — 1) -dimensional boundary of the region B u (r). This is illustrated in 
Fig. [3j The distinction I am making between "hypervolume" and "hyperarea" is somewhat 
arbitrary, since these are both terms involving regions in a higher dimensional space. When 
I wish to talk about a i/-dimensional subregion of the hyperspace M. u , such as B v {r), I will 
use the term hypervolume. On the other hand, when I wish to emphasize a boundary region 
of a hypervolume, such as S v -i(r), I will use the term "hyperarea." Regarding the usage of 
the term "solid angle," suppose we keep the radius r fixed but vary the angles 6$ over ranges 
d9i. The region swept out by this procedure lies on the (y — l)-dimensional sphere S v —\{t) 
with a hyperarea dS u _i = dQ u _i r v ~ x . We are therefore justified in calling dQ u _i the solid 
angle in v dimensions. 

Finally, let us discuss the [y— l)-dimensional cylindrical C u -i(r, L) of radius r and length 
L. Again, it is easiest to argue from analogy in three dimensions. To form a two-cylinder 
C<i{r^ L) in R 3 , we let a two dimensional disk B 2 (r) sweep out a volume as it moves a distance 
L in the orthogonal direction, which is illustrated in Fig. H^i. Similarly, a corresponding 
{y— l)-dimensional cylinder is formed by letting a (y— l)-dimensional ball B v _i(r) sweep 
out a distance L along the orthogonal axis, as illustrated in Fig. [lb. Therefore, the hyperarea 
of the (V — 1) -dimensional cylinder is 



CUi(r,L) = S^ 2 (r)-L. 



(3.19) 



The i/-dimensional hypervolume enclosed by this cylinder is 



V u (r,L) = B v ^{r)-L. 



(3.20) 
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FIG. 5: Definition of the cross section in a general number of dimensions. The incident flux Iq of species a 
is the rate of particles per unit hyperarea normal to the flow. The units of Iq are L 1_I/ ■ T -1 , where L and 
T denote the units of space and time. By definition, the differential cross section da ac is related to the rate 
dR ac , each at angular position Cl, by dR ac (&) — Io da ac (0,). The cross section per unit solid angle about the 
direction is denoted by da a c/dQ. Except for the specification of Cl, this definition does not depend upon 
the dimensionality of space, and the units of da ac are L" _1 . 



B. The Cross Section 



As a physical application of hyperspherical coordinates, let us calculate the form of the 
classical "cross section" in //-dimensions. For simplicity we will consider a projectile striking 
a fixed target, although we can perform a similar analysis in the center-of-mass frame of 
the two particles. Such a scattering experiment is illustrated in Fig. [5], in which a beam 
of incident particles, denoted by the label a, is fired at a target c with incident flux J . 
The rate dR ac (Cl) at which the scattered a-particles enter a given solid angle dQ u -i about 
the direction Q is then measured. The flux Jo is a characterization of the rate at which 
particles move along the beam axis. In v dimensions, the spatial region normal to the axis 
is a [y— l)-dimensional hyperplane, and the flux J is the number of particles per second per 
unit hyperarea passing through this plane. For example, if the beam direction is n, then the 
number of particles in a time interval dt passing through a hyperarea dA^ normal to n is 
given by dN = J ■ dA^ ■ dt. The engineering units of Jo are therefore L x ~ v ■ T -1 . In analogy 
with the usual cross section in three dimensions, we define da ac through 

da ac ■ J = dR ac , (3.21) 

and da ac therefore has engineering units of L v ~ x . 

Suppose the scattering center is a central force, such as the v- dimensional Coulomb po- 
tential. The particle is confined to a two-dimensional plane for central potential motion, and 
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this holds true even in v dimensions. Let b denote the impact parameter of projectile. As 
the particle traverses its plane of motion, its position is uniquely characterized by a function 
b = b(6), where 9 is the angle between the beam direction and the projectile (with the scat- 
tering center defining the origin). From Fig. HJ the number of particles per unit time passing 
through the hyperannulus of width db and radius b is is dN = S u -2{b)db-Io = Vt v ^2 b u ~ 2 db-Io, 
and by particle number conservation, the same number of scattered particles reaches the hy- 
perannulus at 9. The cross section in a //-dimensional central potential is therefore given 
by 

da ac = Q u - 2 b u ~ 2 db . (3.22) 

This is Eq. (8.31) of Ref. 2;], the starting point for the classical calculation. The cross section 
will appear in the Boltzmann equation. To include two-body quantum scattering effects, we 
replace the classical cross section by the quantum cross section: 

i 1 2 d u p c d v p a 
v a — v c da ac = T rr— 7 — rr- i (3.23) 

where T is the quantum scattering amplitude. In the calculations that follow, we shall use 
work in the extreme quantum limit where the Born approximation for the amplitude can be 
employed. 



C. The Coulomb Potential in Arbitrary Dimensions 

Now that we have discussed the cross section in an arbitrary central potential, let us 
concentrate on the special case of the Coulomb potential. The physics of dimensional con- 
tinuation is contained in the //-dependence of the Coulomb potential in //-dimensional space, 
which ensures that short distance physics is emphasized in v > 3 and long distance physics 
in v < 3. Changing the spatial dimension about v — 3 therefore acts as a "physics sieve." 
Let us first construct the electric field of a point charge in v dimensions. Maxwell's equations 
are easily generalized to an arbitrary number of dimensions, and in particular, we can write 

V-E(x) = p(x) , (3.24) 

where E = (Ex,-- - ,E U ) is the electric field vector and V = (d/dxx,--- ,d/dx u ) is the 
//-dimensional spatial gradient. The charge density p has engineering units of charge divided 
length to the u th power, which I will write as Q/L u . In integral form, the equation can be 
written 

J cfxV-E = e , (3.25) 
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where e is the total electric charge contained in the hypervolume E. Note that the dimen- 
sionality of space is now explicitly indicated by the integration measure. We can employ 
the usual symmetry argument to find the electric field of a point source at the origin. Let 
B r be the //-dimensional ball of radius r centered on the point charge e, and denote the 
[y — l)-dimensional hyperspherical boundary by S r . By symmetry, the field E(x) points ra- 
dially outward with a magnitude E(r) along the direction x normal to S r . The length E(r) 
depends only upon the radial distance r = |x| and not upon its angular location along S r . 
The divergence theorem holds in an arbitrary number of dimensions, and since the hyperarea 
of S r is given by f)3.18p . we find: 

e= [ d"xV-E = I dA-E = Q u _ 1 r u ~ 1 ■ E(r) with Q u _ x = 2% - . (3.26) 
J B r Jsr r <>/2) 

The electric field is therefore given by 

E(x) = ^^-x, (3.27) 

where we are using the notation x = r x, with x being a unit vector pointing in the direction 
of x. 

I find it more convenient to work with the electric potential, a scalar quantity 0(r) defined 
by E(r) = —d<ft(r)/dr. In fact, I will work with the potential energy V v = e<f)(r), so that 

, . r(W2-l) e 2 
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where I have appended a subscript to the potential energy to remind us that we are working 
in v dimensions. For two charges e a and e& separated by a distance r, one only need replace 
e 2 by the product e a e&. For u = 3, the geometric factor in (I3.28j) becomes 1/Att, which is the 
origin of the An of rationalized units. Figure [6] shows the Coulomb potential for u = 3, along 
with two representative dimensions on either side of v = 3. As the figure illustrates, the short 
distance behavior of the Coulomb potential becomes more pronounced in higher dimensions, 
while long-distances are emphasized in lower dimensions. For aesthetic reasons, the arbitrary 
integration constant for the potential energy has been adjusted in each case so that all three 
graphs intersect at a single point. This figure illustrates quite dramatically that by simply 
dialing the dimension u, we can dial a potential V v (r) that filters either long-distance or 
short- distance physics. 

In the Born approximation to quantum Coulomb scattering, which we shall employ 
shortly, we need the Fourier transform of the Coulomb potential ( 13. 28ft . Unlike the spa- 
tial representation V u (pc), the Fourier representation of the //-dimensional Coulomb potential 
takes the same form in any dimension, namely, 

K(k) = ~ , (3.29) 

where k 2 = Yli=i * s J us t the square of the norm of the //-dimensional wave number k, 
and I am using the conventions 

K(x) = /^e— k K(k) (3.30) 

K(k) = y^xe ixk K(x) . (3.31) 

With these conventions, the amplitude in the Born approximation in any dimension is given 
by 

e 2 

T B — h— , (3.32) 

where q = p a — is the momentum transfer during the collision. This is a function only 
of the square of its argument q 2 . In particular, the Born approximation does not introduce 
dependence upon the center-of-momentum energy W, and this is what renders its use so 
convenient. 

Expression f!3.29j) for the Fourier transform of the potential (13.281) can be established 
in a number of ways, perhaps the easiest being an straightforward application of Laplace's 
equation, 

V 2 K(x) = e 2 5 M (x) • (3.33) 
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Upon inserting (13.301) for V^(x) into (I3.33j) and using the integral representation of the delta- 
function, we can write Laplace's equation in the form 



d v k 
(2nY 



e~ ix - k k 2 V u {k) =e 2 



d v k 

(27T)" 



-txk 



(3.34) 



and solving for V"(k) provides (I3.29[) . It might also be informative to prove f!3.29|) using 
the more direct approach of performing the Fourier transform directly. Substituting the 
Coulomb potential ( 13.281) into ( 13.311) . and then using ( 13.161) to rewrite the i/-dimensional 
integral, we find 



V u (k) = n v - 2 J drr v - 1 / <lf) sm"- 2 (h 
e 2 r(i//2-l 



v— 2 n Ark cos 9 



r>/2-l) e 2 



2^ r(i//2 - 1/2) 



4^/2 ,. 
rdr I du{l-u 2 )^' 2 



u-2 



^ikru _|_ g— ikru 



(3.35) 



(3.36) 



where we have made the change of variables u = cos 9. It is convenient to keep the exponential 
terms in square brackets rather than converting their sum into a cosine term. We will perform 
the r-integration by deforming the contour slightly off the real axis, 



dr r 



gikru _|_ g— ikru 



co+ie 



drre ikru + 



dr r e 



-i k i 



(ku) 2 ' 

Upon substituting this back into ( 13.361) and changing variables to t = u 2 we can write 



(3.37) 
(3.38) 



K(k) 



e 2 1 I>/2 - 1) 



2^ k 2 T(u/2- 1/2) 



^(i-t)^- 3 )/ 2 r 3 / 2 



(3.39) 



where the second term in the integrand introduces the pole \j[y — 3) into physical quantities, 
and the t-integral takes the form of the Euler Beta function 



B(x,y) 



dtt x ~ x (1 -t) 



2/-1 



r(x)F(y) 



(3.40) 



T(x + y) 

with x = -1/2 and y = u/2 - 1. Using T(-l/2) = -2^ gives (I3T29D . In Section HVB] we 
will need yet another representation of the Beta function, which I record here for convenience: 

T(x)T(y) 



poo 

B(x,y) = / dtt x ~ x {l +ty 
Jo 



-x-y 



T(x + y) 



(3.41) 
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D. Kinetic Equations in Arbitrary Dimensions 



1. Distribution Functions 

A particle in a space M? of arbitrary dimension v G Z + is fully characterized by its position 
and momentum x and p, which have rectilinear coordinates xi and pe for I = 1, ■ • • , v. I will 
often denote the square and the magnitude of the momentum by p 2 = p • p = YTt=\ Pe anc ^ 
p = |p|, respectively. For example, p~ 3 is shorthand for |p|~ 3 = ( Yli=xPe ) ■ A swarm of 
particles distributed over position and momentum values is characterized by a distribution 
function / defined by 

d v x d v p 

/(x, p, t) - — — — = number of particles in a hypervolume d u x about x 

^ 71 ' and d v p about p at time t . (3.42) 

The factor of {2-kK) v in the denominator is a conventional normalization factor, and for a 
spatially uniform distribution f a this gives the normalization 

^UPa) = n a , (3.43) 

where n a is the number density of a-type particles. That is to say, n a d u x is the number of 
particles of species a in a hypervolume d u x, and the engineering units of n a are therefore 
L~ u . From (13.431) . we see that a normalized Maxwell-Boltzmann distribution at temperature 
T a and number density of n a is given by 

/.(P.) =n a ( 2 -^Y /2 ewl-Pa^-} =n a \:e-^ , (3.44) 

where E a = p 2 a /2m a is the kinetic energy and j3 a = 1/T a is the inverse temperature in energy 
units. The thermal wave length for species a is defined by 

X^^m 1 ' 2 . (3.45) 
Consequently, the spatial density of the kinetic energy of species a is given by 

£ >=J^w£r a Mp " th (3 - 46) 

where f a is the corresponding distribution function. 

Suppose now that the rate of change in the distribution function f a is determined by some 
kinetic equation 

^ + v a -V x / a = Y, b K M , (3-47) 
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where V x is the i/-dimensional gradient in position space, v a = p a /m a is the particle velocity, 
and K ab is a scattering kernel between particles of type a and type b. When the distribution 
is spatially uniform we may set the convective term to zero, v a - V x / a = 0, in which case the 
time rate of change in the kinetic-energy density (I3.46P of the a-species is given by 



d£ a f d u p a p 2 a df a ^ f d v p a p 2 a 



We can therefore identify the rate of change in the kinetic-energy density of species a, 
resulting from its Coulomb interactions with species b, by the expression 

' !£a '' f ^ Vl K ab [f] . (3.49) 



dt J (2irh) u 2m a 

Since we are taking each species to be in thermal equilibrium with itself, but not necessarily 
with the other species, each species b is characterized by a unique temperature T b . The rate 
dS a b/dt is therefore proportional to the temperature difference between these species, and 
we write 

^ = -C ab (T a - T b ) , (3.50) 

where C ab is called the rate coefficient. By performing the integrals in (I3.49j) exactly, and 
then comparing with ( 13.501) . we may calculate the coefficients C ab exactly. If the ions are at a 
common temperature T t , then it is more convenient to calculate the rate coefficient between 
the electrons and the collective set of ions, the coefficient C el = J2i C e i of 



2. The Boltzmann Equation 



The derivation of the Boltzmann equation presented in Section 3.3 of Ref. [lOj goes 
through unscathed in a general number of dimensions, and the scattering kernel is completely 
finite when v > 3. The derivation breaks down in v < 3 (that is to say, for v = 1,2,3), 
because in these dimensions the scattering kernel for the Coulomb interaction diverges. This 
behavior for the Boltzmann equation arises because the Coulomb interaction emphasizes the 
short distance physics when v > 3, while the scattering kernel of the Boltzmann equation 
is designed to capture such short distance physics. I will write the Boltzmann equation in 
schematic form as 

^ + v - V x f a = ^2B ab [f} : v > 3 , (3.51) 
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or in explicit form by writing the full scattering kernel as 

Bab[f] = J pjfyv ^ i Vfe ~ V °l {/«(Pa)/ft(Pt) ~ fa(Pa)fb(Pb) j 

(2tt^)-^( P ; + - p a - p 6 ) {2ixh)5(E' a + E' b -E a - E b ) , (3.52) 

with E' c = p' c 2 /2m c and E c = p^/2m c . See Fig. [5] for an explanation of the cross section 
da ah m v spatial dimensions. We can include the quantum effects of two-body scattering, 
to the order in g to which we are working, by replacing the classical cross section by the 
corresponding quantum cross section defined in (13.231) . It then becomes necessary to calculate 
the quantum transition amplitude T(ab — > a'b') = T a / b i. ab , and rewriting (I3.23|) in the form 

|v 6 - v.| - \T a , b , ab \ , (3.53) 

one can then include quantum effects by using the scattering kernel 

Bab[f] = J (2^^ (i^o 17 (l^ 17 1^' 6 '^ 6 ! 2 - 

(2nhy 5"(p' a + p' b - p a - p 6 ) (2nh)S + E' b -E a - E b ) . (3.54) 

For simplicity, in Sec. IIV Al we shall use the Born approximation (13.321) for the transition 
amplitude, which corresponds to taking the extreme quantum limit. When v > 3, expression 
(I3.49P allows us to write the rate of change of the energy density resulting from the now finite 
Boltzmann equation as 

I have used a "greater than" superscript to remind us that we should calculate (I3.55j) in 
dimensions greater than three. 



3. The Lenard-Balescu Equation 

In dimensions less than three one finds a complementary situation to the Boltzmann 
equation, namely, the derivation of the Lenard-Balescu equation is rigorous and completely 
finite when v < 3. This is because the long distance physics of the Coulomb potential is 
dominant in dimensions v < 3, and the Lenard-Balescu equation is designed to capture such 
long distance physics. I will write the Lenard-Balescu equation in schematic form as 

^ + v«- V x f a = J2 b L M -v < 3 , (3-56) 
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where the kernel of the //-dimensional Lenard-Balescu equation is the obvious generalization 
from three dimensions, 

2 



Lab\f] 



d u p a d u p b d v k 



(2-Khy (2-nhY (2tt) 1 



V Pa -k 



e a e b 



fc 2 e(k,v a -k) 



k-V r 



ir8(v a -k - v 6 -k) 



k-V Pa \f a (Pa)f b (Pb) , (3.57) 



where 'V Pc is the //-dimensional momentum gradient. Reference shows that the dielectric 
function of a weakly to moderately coupled plasma is given by 



d u p c 



1 



k ■ V Pc / c (p c ) , 



(3.58) 



k 2 J (2ttH) u uj-k-v c + ir] 

where the prescription 77 — ^ + is implicit and defines the correct retarded response. We can 
use ( I3.58P in (13.571) to the order in g to which we are working. The sum in (I3.58j) is over all 
plasma components, and the velocity v c = p c /m c appearing in the denominator is really an 
an integration variable. Therefore, when v < 3, the rate (I3.49P allows us to express 

d v p p 2 



Lab[f] 



v < 3 . 



(3.59) 



dt J (2nh) u 2m a 

I have used a "less than" superscript to remind us that we should calculate (13.591) in dimen- 
sions less than three. 

It is convenient to express the dielectric function in terms of a complex function F{v) 
defined by the relation 

k 2 e(k, k- v) = k 2 + F(v cos 9) , (3.60) 

where 9 is the angle between k and v. The engineering unit of the argument of F is velocity, 
while the unit of F itself is wave-number squared. Expressions f!3.58j) and (I3.60P imply the 
dispersion relation 

Ptotal(^) 



F(u) 



+00 



dv ■ 

oo v-u-irj 



where the limit 77 — > + is understood, with the spectral weight being defined by 

AotalO) = Y] Pc{v) 

< * r. 



Pc{V) 



2tt 



1/2 



For future use, we shall require the convenient relations 

F(-v) 
lmF(v) 



F*(v 
1 



2i 



F(v)-F*(v) =np tot&l {v) 



(3.61) 

(3.62) 
(3.63) 

(3.64) 
(3.65) 
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While we shall not require the real part 
for completeness I shall record it here: 

ReF(v) = l ~ 

b 

where the Dawson integral is defined by 

daw(x) 

E. Calculating the Rate 

Returning to three dimensions for a moment, and dropping the species index on (13.461) 
for simplicity, the rate of change in the kinetic energy density is simply given by 

f = Jj^L% (p ' t} (3 - 68) 

The problem with a straightforward evaluation of f!3.68|) in three dimensions is that any 
potentially tractable kinetic equation gives a logarithmically divergent result for the rate, 
either at short or long distances, depending on the deficiencies of the particular kinetic 
equation in hand. In principle, however, calculating the rate is a well defined procedure: 
it is the approximation scheme employed in finding the requisite distribution function that 
introduces a divergence. In other words, if one knew the exact single-particle distribution 
function f(p,t), then the rate (13.681) would be finite. However, the one-point function 
fi = f(Pi t) of the BBGKY hierarchy can be known exactly only by solving the entire set 
of coupled multi-particle correlation functions exactly, an impossible feat. Hence, we must 
approximate the exact distribution function by one obtained through truncating the BBGKY 
kinetic equations. Even worse, the truncation process is rather subjective in that it depends 
upon the type of physics one deems important; for example, truncation to the Boltzmann 
equation is only useful if we can neglect long distance correlations (which, in this problem, we 
cannot 9 ). Conversely, truncation to the Lenard-Balescu equation captures the long distance 
physics, but misses the short distance physics. In a nutshell, then, our problem is the 
following: to calculate the rate, we require the exact distribution function of the full hierarchy 
of kinetic equations, a problem we cannot hope to solve without a Quantum Computer or 

9 In reducing BBGKY to the Boltzmann equation, we make the approximation that two-body collisions 
are uncorrelated, thereby allowing the replacement fi — > fi ■ fi in the scattering kernel. However, the 
correlations described by /2 and higher act back on fi to render the integrals in (|3.68p finite at long 
distances (and such correlations are neglected by the Boltzmann equation). 



of F, nor is there space to compute this function, 




(3.66) 



(3.67) 
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a Mentat; we must therefore truncate the kinetic equations, but this does violence to either 
the short or long distance physics, thereby introducing spurious divergences. 

A hint out of this Catch 22 comes from the following observation. The truncation problem 
only occurs for the Coulomb potential, and only then in three spatial dimensions. As we 
saw in (13.281) . the Coulomb potential in v dimensions is V v {r) = C u e a eb/r u ~ 2 , where C u = 
T(y/2 — 1)/Ati u / 2 is a geometric constant, and this form of the potential renders the scattering 
kernels of the Boltzmann and Lenard-Balescu equations finite, except for the single case v = 3 
(ironically, the case of interest). Let us therefore start with the exact BBGKY hierarchy in 
v spatial dimensions, with the understanding that taking v ^ 3 is a regulating procedure, 
and that we must eventually return to three dimensions. This procedure, however, is robust 
enough to capture the correct physics as the limit v — > 3 is taken. 



1. Reduction of BBGKY 

In the exact same manner as in the last few section, one can generalize the BBGKY hier- 
archy to an arbitrary number of dimensions. Furthermore, as we discussed in Lecture I 0], 
when the number of spatial dimension is greater than three, BBGKY reduces to the Boltz- 
mann equation (13.511) and (13.521) to leading order in the plasma coupling g. Conversely, 
when v < 3 the BBGKY hierarchy reduces to the Lenard-Balescu equation (13.561) and (13.571) 
to leading order in g. As discussed in the previous section, these reduced kinetic equations 
(Boltzmann and Lenard-Balescu) are finite in their respective dimensional regimes. In other 
words, besides rendering the truncation process finite, the physical utility of keeping the 
dimension of space arbitrary is that in dimensions greater than three, the leading order in g 
behavior of BBGKY is just the finite z/-dimensional Boltzmann equation: 

BBGKY in u > 3 ^ + v - V x f a = J2 b B M to LO in g , (3.69) 

where B ab [f] is given by (13.521) . As discussed at length in Lecture I, expression (13.691) is the 
point at which the physics of dimensional continuation enters the calculation: dimensions 
greater than three select for short distance physics. For v > 3, the rate of energy transport 
from plasma species a to species b is therefore given by the finite expression (13.551) . Turning 
now to dimensions less than three, we have seen in Lecture I that the leading order in g 
behavior of BBGKY reduces to the Lenard-Balescu equation, 

BBGKY in v < 3 ^ ^ + v a -VJ a = J^M/] to LO in g , (3.70) 

where the finite i/-dimensional scattering kernel L a b[f] is given by (I3.57p . Again, the physical 
content of dimensional continuation enters at this stage: dimensions less than three select 
for long distance physics. This dimensional reduction of BBGKY is illustrated in Fig. [3 
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FIG. 7: For v>3 the "textbook derivation" of the Boltzmann equation for a Coulomb potential is 
rigorous; furthermore, the BBGKY hierarchy reduces to the Boltzmann equation to leading order 
in g. A similar reduction from the BBGKY hierarchy holds for the Lenard-Balescu equation in 
v < 3, and the "textbook derivation" is also rigorous in these dimensions. In v = 3, the derivations 
of the Boltzmann and Lenard-Balescu equations break down for the Coulomb potential. 

From the rate equations (I3.47p - fl3.49p . in dimensions v > 3 expression f!3.69j) gives the 
rate of energy transport from species a to species b as 

dpt = f <?Va Pl_ b [f] LQ in (3 71) 

dt J [2ttK) u 2m a UJ y ' v ' 

which, as we shall see, takes the form 

dS > a 2 

— ^ = H{y) + 0(u - 3) to LO in g when v > 3 . (3.72) 

at v — 3 

We have omitted the species indices from H{y) for simplicity, and we shall calculate this 
quantity in Section TlV A[ Rather than a logarithmically divergent result, we obtain a finite 
answer involving a simple pole of the form \j[y — 3), which of course is the origin of the 
divergence in three dimensions. Since we will eventually return to v = 3, there is no need to 
calculate the 0{v — 3) terms in f)3.72p . as these terms vanish when v — > 3. Similarly, from 
f)3.70p the corresponding rate in energy transport from species a to b is therefore 



toL0 " 19 ' (373) 



dt 

From the calculation in Sec. IIVBI we shall find 



dS < a u ~ l 

-J* = G{v) | + 0(3 - v) to LO in g when v < 3 , (3.74) 

ox o jy 
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where we have omitted the species indices from G(u) to save writing. Note that the leading 
behavior of both (13.721) and (13.741) is formally of order g 2 in three dimensions, illustrating that 
neither short nor long distance physics dominates in v = 3, but rather, that ultraviolet and 
infrared length scales contribute to the same order in three dimensions. It is a property of the 
Coulomb potential itself that v = 3 is the fulcrum around which the short and long distance 
physics pivot. Using the appropriate Coulomb potential (13.281) for V v {r) in the scattering 
kernels of ( 13.711) and ( 13.731) . the integrals now converge, and they are calculated exactly in 
Sections 7 and 8 of BPS [2]. In these notes we shall calculate them in Sections IIV Al and II V Bl 

2. Obtaining Next-to- Leading Order from Leading Order 

In three dimensions, or in the limit v — > 3, we still are still plagued by the long and short 
distance divergences from the simple poles in ( 13.721) and ( 13.741) . a problem we must now 
confront if we are to obtain a meaningful result. As described in Lecture I, to compare the 
rates (I3.72p and (I3.74p . we must analytically continue one or the other to a common value of 
the dimension v. Analytically continuing the spatial dimension makes sense because we can 
view the quantities d£> b /dt and d£< b /dt as functions of a complex parameter z/, even though 
they were only calculated for positive integer values of v. This is analogous to analytically 
continuing the factorial function on the positive integers to the Gamma function on the 
complex plane. For defmiteness, I will analytically continue (I3.74p to v > 3, in which case 
the ^-dependence becomes subleading relative to the g 2 dependence of (I3.72p . The analytic 
continuation of ( 13.741) takes the same functional form for any v 6 C, but in this section I 
will write the analytic continuation as 

d£< 2+(^-3) 

= -G(u) - — + 0(u - 3) to NLO in g when u > 3 . (3.75) 

(Xv is o 

Since we are now working in the regime v — 3 > 0, I have written the exponent of g in a 
form to emphasize that g 2 ^> g v ~ l when g <C 1 and v > 3. There are no terms with powers 
of g intermediate g 2 and g"' 1 , so the analytic continuation of dS^ b /dt to parameters v > 3 
is not only subleading in g, but it is indeed next-to- leading order relative to ( 13. 72ft . This is 
illustrated in Fig. [HJ 

In the rates (I3.72p and (13.741) . we need to work consistently only to linear order in the 
small parameter e = v — 3; therefore, we should expand H{y) and G{y) to first order in e, 
allowing us to write 

H(v) = -A + eHi + Oie 2 ) (3.76) 
G{u) = -A + eGi + 0(e 2 ) . (3.77) 
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analytically continue 
around the v = 3 pole 




d^/dt d£<Jdt 



g v-l = g 2-(3-u) 3 gV -\ _ g 2+(v-3) 

„2-|f-3| „2+|i/-3| 

LO: large when g <C 1 NLO: small when g -C 1 



FIG. 8: The analytic continuation of d£< b /dt from ^ < 3 to the region > 3 in the complex j/-plane. The 
same expression can be used for d£^ b /dt throughout the complex plane since the pole at v = 3 can easily 
be avoided. The quantity d£^ b /dt ~ y 2 +( ,y - 3 ) is leading order in g for v < 3. However, upon analytically 
continuing to > 3 we find that d£^ b /dt ~ g 2 +l"~ 3 l is next-to-leading order in g relative to d£^ b /dt ~ g 2 . 

It is crucially important here that and G(u) give the same value at v = 3, a term 

that I have called A in (13. 76ft and ( 13.771) . otherwise the divergent poles will not cancel. In 
Section HVl we shall calculate the coefficients H(u) and G{y), and we will indeed explicitly 
see that H(y = 3) and G[y = 3) are equal. We will also calculate H\ = H'{y = 3) and 
G\ = G'(u = 3) in closed form, thereby providing an exact result for the rate to leading and 
next-to- leading order in g. 

Finally, upon writing g e = exp{elng} in (I3.75p . and expanding to first order in e, we find 

^ = - + \ng + 0(e) . (3.78) 
e e 

This is where the nonanalyticity in g arises, i.e. the \ng term, and we can now express the 
rates (13TT2D and (13TT5D as 

df > A 

g 2 + H ig 2 + 0(u-3;g 3 ) u>3 (3.79) 



dt v — 3 

dS< b _ A 



dt v — 3 



g 2 -G ig 2 -Ag 2 \ng + 0(u-3-g 3 ) v>3 . (3.80) 



These expressions hold in the common dimension v > 3, and to obtain the leading and 
next-to-leading order result in three dimensions, we add and take the limit: 



— - — = hm 

dt y^3+ 



dt dt 



0(g 3 ) = -Ag 2 Ing + Bg 2 + 0{g 3 ) , (3.81) 



with B = H\ — G\. Compare this with the rate ( 12. 21) . or the alternative expression ( 12. 31) . 
In this way, BPS calculated the energy exchange rate from Coulomb interactions between 
plasma species, accurate to leading order and next-to- leading order in g. 
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IV. CALCULATING THE RATE TO SUBLEADING ORDER 



We now perform the real calculation of the energy exchange rate between electrons and 
ions. BPS |2J considered the general case, finding C a b in (I3.50P for any collection of plasma 
species in any quantum regime. In these notes, however, I will only consider the extreme 
quantum limit valid at high temperatures. This is the case of most general interest, and it 
is also the case in which the algebra simplifies considerably. We shall also take the electrons 
to be in equilibrium with themselves at temperature T e and the ions in equilibrium with 
themselves at temperature T r , another situation of general interest. Upon summing over the 
ions, the rate equation then becomes 

d£ e i 



dt 

where the collective rate coefficient that we shall calculate is 



-C el (T e -T x ), (4.1) 

we shall calculate is 

C eI = V.C ej . (4.2) 
< '1 



A. Boltzmann-Equation: Short-Distance Physics 

We now work in v > 3 dimensions where the short- distance physics dominates. To 
calculate the rate of change of the electron distribution, we will employ the Boltzmann 
equation with two-body quantum effects in the Born approximation. The rate of energy 
exchange from the electrons to the ions is therefore, 

ggg / d v p e yj dfejpe) f d v p e p\ 

dt J (2nhY2m e dt (2-KhY 2m e ei[Jl ' 1 1 

where we have taken a = e and b = i in (13.511) and (13.541) . and the electron and ion 
distribution functions are given by (13.441) . Using the crossing symmetries p e <-> p' e and 
Pi Pi °f t ne scattering amplitude T B in (I3.54p . the rate of energy exchange from the 
electrons to an ion species i can be written 

d£> _ f d»p' e d% d»p e d v Pi , |2 p' e 2 - p\ f ,. ff . 

J B Je{Pe)Ji{Pi) 



dt J (2nh) u (2nhY (2nhy {2-kK) v 11 2m e 

(W <P5 + P'e - P. - P«) M S + ^f) , (4.4) 

where we dropped the gradient term in (I3.54p because of spatial uniformity. Summing over 
all ions in (14.41) gives the total rate (14.31) . We define the momentum transfer q and the 
average momentum p of the initial and final electron momentum, 

q = Pe - Pe = Pt ~ P t (4.5) 
1 



P 



P ; + Pe . (4.6) 
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Upon performing the ^--integral using the momentum conserving delta-function to set 

Pi = Pi + Pe - Pe = Pi - q , (4-7) 



and expressing the electron momenta in terms of q and p 



1-1 , . 

Pe = P+ 2 C l Pe = p-~Cl, (4.8) 

we can simplify ( 14.41) to read 

<9£> f d v p d v q \ . .\imi _ 

/(^(^(^-S^-^ 2 )^- (4 ' 9) 

We have used the fact that the energy conserving delta function and the energy loss factor 
become, 

ti^ti = — p-q (4.10) 

2m e m e 

5 (Pl^A + Pl^A) = 5 (H-m + 7 pl . (4.11) 

\ 2m e 2mj / \ m e 2mj 2rrii J 

We now perform the ^-integration. Since we will find a similar integral in the next 
section, I will perform a more general calculation here. There will be times when we need 
to integrate a Gaussian and a delta-function, which I will write as 

d Pb r,t_ Tfl ,-O.Eh. to I 1 



-5{k-v b -V)e~^ Eb = hp b (^) e -5^ y2 = ^A 1 ^e-^" m ^ 2 ,(4.12) 



where k is a fixed unit vector (typically another integration variable), and V is a scalar 
independent of p b = rribVb. The integral (I4.12p will be required in several places throughout 
the text, so we will perform the calculation here. Since V = v& • k defines the component 
of V;, parallel to k, will decompose the integration variables vj, into parallel and normal 
components 

v 6 = v ± + (v 6 -k)k = v ± + Vk . (4.13) 
Since v\ = v\ + V 2 , we can write 

= (|H) Y rV«^ W ' f°° dV<S(V - V) e -'it , * m > v * (4.14) 
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and this yields expression (I4.12p . As an application, we will frequently run across an integral 



d u Pb 



6(v b -k-V)f b (p b 



n b 



f (3 b m b \ 

V 27T J 



1/2 



e -\(3 b m b V 2 



(4.16) 



which follows directly from (I4.12p . Using the integral (I4.16p . and taking f!3.44j) for / e (p — q/2) 
gives 



dt 



m P 



(27r/3 i m i ) 1/2 A 



d u p d u q 



-n|T B (g)| 2 p-q 



exp 



Pim 
2ml 



m P 



mi 



Pe 



(p-q) + — 1 - — gpq 



Pi 



Pe 



2m P 



p z + 1 + 



m e Pi\q 



rriiPe I 4 



(4.17) 



When P e = Pi, the linear term in the exponential involving p ■ q vanishes. The integrand is 
even in both p and q, except for the prefactor p ■ q; therefore, keeping q fixed and integrating 
over p gives zero, 



d£> 



dt 







as it should for equal electron and ion temperatures. 

We now examine the general case when the electron and ion temperatures differ. Com- 
pleting the square for the terms in the first square-brackets suggests changing variables to 



P 



m P 



rrii 



Pe\ q 



P+— l-ir i 



Pi 



(4.18) 



and dropping g 2 -terms that are down by relative factors of m e /rrii gives 



dt 



n P rii 



m P 



(27ip tmi ) 1/2 X u 



d u p' d v q 



exp< — 



(2nh) v (2irh) 

Pirn , 2 P, 
(p -q) 



rrii \Pi 



^+(p'.q) 2 + | 
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(4.19) 



2m\ vx ~ ±y 2m e 

where we have expanded the new integration variable p' into normal and parallel components 

p' = p' ± + p' q with f = p q . (4.20) 

Finally, we note that the term (p'-q) 2 in the square brackets of the exponential is down by 
a factor m e /mj relative to the first such term, and we can write (I4.19P as 



dt 



n e rij 



(2np imi ) 1/2 XI 



dF 



rrii \P: 



2nh J (2nh) v - 1 (2nh) 
exp 



Pi m i p/2 Pe 



2ml 



2m P 



-,2 , T 
TP- +T 



(4.21) 
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The linear P' term in the prefactor integrates to zero, and we arrive at 



dt 



2m,- 



(2nf3 t m t y /z \ 



1/2 ,„ I Pe 



Pi 



dF f <T- l p' ± d»q ,,2 

ft T b(<?) 



2nh J {2nhy- 1 (2nh) u 



i _ fiiTflj p/ 2 _ Pe 

2m 2 2m e 



-,2 . 9 

P± + ~ t 
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The integrand is now Gaussian in the variables P' and p{ 2 , and we find 



(27r/3 i m l ) 1/2 

K 

Substituting back into (14.221) gives 
dS> { n e rii X e m e ( (3 e 



OO JTTD/ 



dF' 
2nE 



cxp 



2ml 



P'^ = — 

h 



d v - l p' 
(2nhy 



- exp 



Pe -i2 



2m, 



dt 2m,- h 



l 



d u q 
(2irH)' 



h \T B (q)\ q exp< - 



Pe q 2 



2m e 4 



(4.22) 

(4.23) 
(4.24) 

(4.25) 



Up to this point, the two-body scattering amplitude T could have been general, but we now 
explicitly employ the Born approximation 

eei 



so that 



d v q 1 



{2ithy q-- 



exp< — 



Pe q 2 



2m P 4 



(4.26) 



(4.27) 



where the ion plasma frequency is uof = ejrii/mi and the electron Debye wave- number is 
Kg = e 2 n e /T e . 

Another expression we will encounter is the gamma-function, 



T(z) = / duu z -'e~ u Re(» > 



(4.28) 



o 

The gamma-function (I4.28P provides a nice trick for calculating integrals of the form 

"OO 



a 



n—l „—a i 



ds s n ~ l e 



(4.29) 



r(n) Jo 

In doing calculations, we will often need to take the norm of a vector to some power, such 
as |q|~ < ^~ 3 ' ) / 2 . For example, we can use (14.291) to exponentiate the norm into a more easily 
handled Gaussian by writing |q| = (q • q) 1 / 2 , so that 

, „ 1 



T(m/2) 



dss 771 ' 2 - 1 e~ sqq . 



(4.30) 
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We can now perform the g-integral to give 



dV( l g -3 c -/W/8m e 



frcWJo dqq 



fie 2 
S = a ? 



(4.31) 



2tt^ 2 



(27rft)" r(i//2) V A 



5771, 



(v-3)/2 



1 T(3/2) / 4 \(^ 3 )/ 2 ^ /„_ 3 



dss ( " _5)/2 e- s (4.32) 



r 



47r 2 n 3 I>/2) VA 2 y " V 2 

where T(3/2) = a/vt/2, and substituting this back into 1 14. 27ft allows us to express 



d£> k\ (\ e m e x 2 



^ 2tt V 27T/7 



\ T(3/2) 1 / 4 



JT(u/2) 2 Va 



("-3)/2 



i/-3 



Note that 



Km. 



(3 e m t 



1/2 



(4.33) 



(4.34) 



(4.35) 



27T/7 V 27r 

and therefore the rate coefficient becomes, upon dividing by the temperature difference and 
then summing over the ion components, 



°el 



K 



fi e m e \ 1/2 T(3/2) 1 (A 



2tt 1 V 2tt J T(u/2) 2 \A 



("-3)/2 



r 



i/-3 



(4.36) 



B. Lenard-Balescu Equation: Long-Distance Physics 



We now wish to calculate the leading order long-distance physics by working in spatial 
dimensions v < 3. The rate of energy exchange from the electrons to the ions is 

d£e = 2 f dV Pe Vl dfJjPe) , , 

dt J (2-KhY 2m e dt ' 1 ' 1 

and from 113.571) we find the rate of energy exchange between a spatially uniform distribution 
of electrons and ion species i, 



d lA = _ 2 f dUp * dUp * duk P * V „ k 



(2nh) u {2-nhY (2tt)^ 2m e Pe 



7T<5(k-v e — k- Vj 



k 2 e(k, k- v e 

k- Vi A-k-VeAe)/ e (Pe)/i(Pi) ,(4-38) 



where we have used the distribution (13.440 to write 'V Pb fb(Pb) = — fib v & /b(Pfc)- Integrating 
p e by parts using k • V Pe (p 2 /2m e ) = k • v e , and integrating over the ion distribution with 
f l4TT6|) gives 

d u p e d u k 7r e 2 k-v e pj(k • v e 



at 



(27T/z)- (27r)" |A;2 + F(k-v e 



/e(Pe) , 



(4.39) 
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where we have used (13.601) and (13.631) . When the electron and ion temperatures are equal, 
(3 e = Pi, we see that the rate vanishes, as it must. Inserting unity in the form 

/CO 
dv6(v -k-v e ) = 1 (4.40) 
-oo 

allows us to express the rate as 

~(Pe «\ f°° , f d u p e d u k ire 2 kvpi(v) p , * , , 

"af = 2 (s ~ l )L dv J^Wr\i^ s(v - k - v ' ),M • (441) 

and upon performing the electron momentum integrals with ( 14. 16ft . we find 



9* V 2tt ; I V J (2tt)- /_„ |*HF(b)| j 



where we have used (11.31) for K e . 

For individual ion species i, the ^-integral can only be performed numerically; however, 
if we sum over all ion species we can perform the integral by completing a contour in the 
complex t>-plane JsJ: 

/OO / \ 2 

= ,4,3) 



Summing over the ions in (14.421) therefore allows us to express 



2K>f (^) I/2 ( Tl _ Te )./^_±_. (4 . 44) 



Of 

Note that the fc-integral 

(27r)" (P + K 2)2 (4 ' 45) 

converges both in the UV and IR. This is because the K 2 e term in the denominator provides 
a long-distance cutoff (small values of k), while the condition v < 3 provides short-distance 
convergence (large values of k). This integral can be performed by first converting to hyper- 
spherical coordinates and then changing variables by k 2 = k 2 t: 



e 
oo 



; k^dk — k ^ =^-^± / dtt-^it + iy 2 , (4.46) 



(27r)-y (k 2 + K 2 e y 2 (2tt) 

where the solid angle integrals Q u -± are given in ( I3.26p . The t- integral can be recognized as 
the Beta function, defined in (I3.4ip with x = (y + l)/2 and y — (3 — f)/2. Inserting the 
appropriate Gamma functions and factors of 7r gives 

« t _ i rp/2)/^ fiV-* , (4 . 47) 



(2tt)^ (P + K 2)2 4^2 r(z//2) V 2 y V 2 J \4tt 
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where T(3/2) = Substituting (1PT1) back (Oil) gives 

or the rate coefficient 

cl ^(f^f^M—U—) ( Kj V )r '- («») 

61 2tt 1 V 2tt y r(z//2) 2V2yv2yv 47r y 

C. Combining the Leading and Next-to-Leading Order Terms 

Recall that the rate coefficient in v < 3 and v > 3 takes the exact analytic form 

c « = s (— ) fm2 T {—) UJ (4 ' 50) 

= * f/^V /2 £2a I r f^ r (i^j . (4 . 51) 

2tt 1 V 2vr y r(i//2) 2V2yV2yv 4 W 

The gamma functions in the above expressions can be expanded in the parameter e = v — 3: 

r(e) = i- 7 + 0(e) (4.52) 
e 

r(l + e) = 1 - 7 e + 0(e 2 ) . (4.53) 

We often need to multiply a term of the form A u ~ 3 by a pole l/(z/ — 3). This will produce 
a pole term and a finite contribution. In fact, this is the origin of the coefficient under the 
logarithm, and since this point is so important, I will reiterate it once again. For ease of 
notation let e = v — 3, so that the v — ► 3 limit is the same as the e — > limit. In any 
calculation we can therefore drop terms 0(e) and higher; however, we must be careful and 
drop such terms too soon. This is because an order e term could multiply a pole term of the 
form 1/e, thereby giving a finite nonzero result in the limit e — > 0. The following example 
illustrates this point. Let us consider the product of A 6 with the pole 1/e. We first expand 
A e in powers of e as follows 

A e = exp{ln A e } = exp{e In A} = 1 + e In A + C(e 2 ) . (4.54) 

Upon multiplying this expression by the pole we find 

= ~ + lnA + 0(e) . (4.55) 
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Therefore, upon using expression (14.551) in (I4.50p and (I4.5ip . we find 



v + 1 



r 



z/-3 
2 

3-v 



This gives the rate coefficient 



Cei — + C 



< 

ei 1 ""el 



2 

— W T 

2tt 1 



" 4 

4 

47T 



2tt 



(f-3)/2 



(u-3)/2 



u-3 
2 

z/-3 



4 

+ A? 



7 



In 



{4^} 



(4.56) 
(4.57) 



1/2 



In 



167T 



-7-1 



(4.5* 



Note that the pole terms have canceled, rendering a finite result accurate to leading and 
next-to- leading order in g. As we have seen, the next-to-leading order term gives the exact 
coefficient under the logarithm (including the term —7 — 1). The argument of the logarithm 
can be expressed as 



16tt 8m P T P T e 



8T 2 



ft 2 e 2 n P h 2 lu 2 



and therefore the rate coefficient takes the form 



(3 e m e 



In A. 



In 



2vr 
8T 2 



h 2 uj 2 



InAp 



-7-1 



(4.59) 

(4.60) 
(4.61) 
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